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8  Different  dislocation  processes  are  shown  to  be  operative  under  high  rate  loading  by  impact-induced 

9  shock  tests  as  compared  with  shockless  isentropic  compression  experiments  (ICEs).  Under  shock 

10  loading,  the  plastic  deformation  rate  dependence  of  the  flow  stress  of  copper  is  attributed  to 

11  dislocation  generation  at  the  propagating  shock  front,  while  in  shockless  ICEs,  the  rate  dependence 

12  is  attributed  to  drag-controlled  mobile  dislocation  movement  from  within  the  originally  resident 

13  dislocation  density.  In  contrast  with  shock  loading,  shockless  isentropic  compression  can  lead  to 

14  flow  stress  levels  approaching  the  theoretical  yield  stress  and  dislocation  velocities  approaching  the 

15  speed  of  sound.  In  iron,  extensive  shock  measurements  reported  for  plate  impact  tests  are  explained 

16  in  terms  of  plasticity-control  via  the  nucleation  of  deformation  twins  at  the  propagating  shock  front. 

17  ©  2009  American  Institute  of  Physics.  [DOI:  10.1063/1.3067764] 

18 


19  I.  INTRODUCTION 

20  In  the  1980s,  Follansbee  et  al.]  reported  a  strong  upturn 

21  in  the  flow  stress  of  copper,  determined  at  0.15  strain  and 

22  300  K,  when  tested  at  the  upper  end  of  higher  strain  rates 

23  achieved  in  split-Hopkinson  pressure  bar  (SHPB)  experi- 

24  ments.  The  resulting  flow  stress  dependence  is  shown  in  Fig. 

25  1.  A  number  of  other  metals  were  observed  to  have  similar 

26  upturns  in  flow  stress  at  strain  rates  of  the  order  of  104  s_1. 

27  Figure  1  also  displays  “asymptotes”  for  the  low  and  high 

28  strain  rate  dependencies  of  the  activation  area  for  a  thermally 

29  activated  dislocation  model  for  plastic  flow.  In  the  figure,  the 

30  asymptotes  are  labeled  with  the  limiting  values  of  the  acti- 

31  vation  area  A*  determined  from  the  expression 


32 


A*=—[ 9  ln  A 

b  \  dr*  )j 


(1) 


33  where  k  is  Boltzmann’s  constant,  T  is  absolute  temperature, 

34  b  =  0.255  nm  is  the  dislocation  Burgers  vector  of  copper,  y  is 

35  the  shear  strain  rate,  and  f*  is  the  thermal  component  of 

36  shear  stress. 

37  For  connection  of  Eq.  (1)  with  the  material  parameters  in 

38  Fig.  1,  At*  =  Acr/m,  e=  yl in,  and  m=  3.08  is  the  Taylor  ori- 

39  entation  factor  relating  single  crystal  shear  stress  and  strain 

40  to  the  corresponding  polycrystal  effective  stress  and  strain; 

41  see,  for  example,  the  thermal  activation  description  for  dis- 

42  location  dynamics  given  by  Armstrong."  In  accordance  with 

43  that  model  description  of  the  polycrystal  flow  stress  cre  at 

44  constant  tensile  or  compressive  strain  e, 

45  cre  =  mr=  m(  f  +  tc  +  ks£~112),  (2) 


46  where  rG  is  an  athermal  shear  stress  that  is  dependent  on  the 

47  dislocation  density  and  solute  concentration,  ks  is  the  so- 

48  called  microstructural  shear  stress  intensity  for  overcoming 
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the  grain  boundary  resistance,  and  €  is  the  polycrystal  grain 
diameter. 

As  will  be  seen  for  iron,  the  slip-determined  kE=mks  and 
a  larger  twinning-determined  kT  play  important  roles  in  de¬ 
termining  the  initial  yielding  responses  to  shock  loading.  The 
evaluation  of  Eq.  (2)  at  F~1/2  =  0  is  taken  to  specify  the  aver¬ 
age  dislocation  friction  stress  cr0=m(T*  +tg).  Also  for  iron, 
the  temperature  and  strain  rate  dependencies  of  crE  are  in  <x0. 
Otherwise,  as  noted  in  Fig.  1  for  copper,  a  relatively  large 
value  of  ,4*~  1000/?2,  corresponding  to  (Acr/Alne)r 
—  0.8  MPa,  is  measured  as  a  typical  value  for  a  face- 
centered-cubic  (fee)  metal  tested  at  conventional  tensile  or 
compressive  strain  rates,  and  such  measurements  are  gener¬ 
ally  attributed  to  the  thermally  activated  overcoming  of  dis¬ 
location  intersections  during  the  slip  process. 

In  the  present  article  we  account  for  the  strong  upturn 
shown  in  Fig.  1  for  the  copper  flow  stress  measurements  by 
making  connection  on  a  dislocation  nucleation  basis  with 
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FIG.  1.  The  flow  stress  dependence  on  strain  rate  of  oxygen-free  electronic 
copper  at  0.15  strain  and  300  K  as  reported  by  Follansbee  et  al.  (Ref.  1)  and 
fitted  with  asymptotic  activation  area  A*  values  at  the  low  and  high  strain 
rates. 
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other  sharply  rising,  shock-induced,  plastic  flow  stresses 

68  measured  both  for  copper  by  Swegle  and  Grady  and  for  iron 

69  by  Arnold.4  The  upturn  in  the  SHPB  measurements  of  copper 

70  by  Follansbee  et  al. 1  is  shown  to  have  provided  early  evi- 

71  dence  pointing  toward  the  substantially  higher  shock- 

72  induced  plastic  flow  stress  measurements  reported  by  Swegle 

"I 

73  and  Grady.  The  result  is  attributed  to  a  changeover  from  the 

74  often  employed  Orowan  equation,  for  relation  of  dyldt  and 

75  the  dislocation  velocity  v  in 


76 


(3) 


77  where  p  is  the  dislocation  density  to  the  alternative  disloca- 

78  tion  nucleation  equation  for  dyldt ,  also  obtained  from 

79  Orowan5  as 


80 


dy 

dt 


(4) 
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FIG.  2.  Flow  stress  dependence  on  strain  rate  reported  by  Follansbee  et  al. 
(Ref.  1)  (triangles)  and  shock-induced  plasticity  results  reported  by  Swegle 
and  Grady  (Ref.  3)  (circles),  compared  with  ZA  thermal  activation  model, 
the  Swegle-Grady  empirical  relation,  and  the  thermal  activation  dislocation 
nucleation  model. 


81  where  dp!  dt  is  the  rate  of  nucleation  and  Axd  is  a  small 

82  distance  associated  with  the  nucleated  dislocation  displace- 

83  ment. 

84  Arnold4  measured  the  Hugoniot  elastic  limit  (HEL)  pres- 

85  sure  and  follow-on  jump  to  a  higher  plastic  flow  pressure 

86  exhibited  for  plate  impact  experiments  conducted  on  differ- 

87  ent  grain-sized  iron  materials.  The  tests  covered  a  wide  range 

88  of  projectile/target  thicknesses.  Comparison  of  the  results 

89  with  postshock  metallographic  observations  and  with  the 

90  previously  mentioned  grain  size  dependencies  for  slip  or  de- 

91  formation  twinning  in  iron  provides  evidence  for  roles 

92  played  by  both  deformation  responses  and  transition  between 

93  them  in  determining  reasonably  strong,  although  more  so  for 

94  twinning  than  slip,  grain  size  dependences  of  the  HEL.  On 

95  the  other  hand,  because  a  substantial  shear  strain  is  imposed 

96  at  all  points  along  a  propagating  shock  front  with  consequent 

97  relaxation  produced  in  iron  by  nanoscale  twin  nucleation,  the 

98  higher  shock-induced  plastic  flow  stress  will  be  shown  to  be 

99  unaffected  by  the  material  grain  size. 

100  Last  but  perhaps  most  interesting  is  an  interpretation  put 

101  forward  for  recent  isentropic  compression  experiment  (ICE) 

102  results  reported  for  copper  by  Jarmakani  et  al.6  No  shock  is 

103  produced  in  the  uniform  loading  of  an  ICE  test,  and,  there  - 

104  fore,  the  flow  stress  levels  and  dislocation  velocities  are 

105  shown  to  be  limited  only  by  the  energy  dissipation  resulting 

106  from  interaction  of  the  moving  dislocations  with  lattice  vi- 

107  brations,  that  is,  the  so-called  influence  of  “phonon  drag.” 

108  Consequently,  the  flow  stress  may  increase  to  the  theoretical 

109  cohesive  limit  with  velocities  approaching  the  speed  of 

110  sound.  Quantitative  evaluation  of  the  results  of  Jarmakani  et 

111  al.6  leads  to  the  conclusion  that,  at  the  highest  peak  pressure 

112  reported  of  —52  GPa,  the  plastic  flow  stress  is  very  near  to 

113  the  theoretical  limit  for  copper. 


Grady.3  At  the  relatively  lower  stresses  shown  for  the  SHPB  12® 
results,  the  data  follow  a  thermally  activated  dislocation  me-  121 
chanics  based  relationship  proposed  by  Zerilli  and  Arm-  122 
strong,  as  updated  recently  by  Zerilli  in  the  fee  form  123 

<tf.  =  BQem  exp(-  (A,  -  A  In  e)T)  +  crc  +  kj~v2  (5)  124 


in  which  the  experimental  constants  crG,  B0,  /30,  ,  and  ke  125 

were  determined  from  other  experimental  results.  The  first  126 
term  on  the  right-hand  side  of  Eq.  (5)  is  a  =  in  t  in  accor-  127 
dance  with  Eq.  (2).  The  empirical  relationship  128 

cr^em  (6)129 

a 

proposed  by  Swegle  and  Grady  is  also  shown  in  the  figure  130 
as  the  dashed  line  passing  through  their  pair  of  shock-  131 
induced  plasticity  measurements.  The  filled-black-circle  132 
stress  values  are  transformed  from  the  measured  shock  pres-  133 
sures  P  using  the  relation  134 


<r  = 


1-2^ 


- P, 

l-v 


(7) 


135 


where  v  is  Poisson’s  ratio.  A  value  of  Poisson’s  ratio  v  136 
=  0.345  for  copper  was  used.  137 

Q 

Armstrong  et  aV  proposed  that  the  boundary  conditions  138 
at  the  propagating  shock  front  required  sequential  generation  139 
of  a  nanoscale  dislocation  structure  so  that  Orowan’s  second  140 
equation  [Eq.  (4)  above]  should  be  the  appropriate  relation-  141 
ship  to  employ  in  a  thermal  activation  model.  Further,  as  142 
seen  from  the  result  shown  in  Fig.  1,  a  lower  limiting  value  143 
of  the  area  of  activation  of  order  of  ~b2  dimension  applies.  144 
By  incorporation  of  a  lower  limiting  activation  volume  V0,  145 
then  in  a  thermally  activated  relationship  for  dpldt,  the  146 
simple  relation  is  obtained  147 


114  II.  CONVENTIONAL  STRAIN  RATE,  SHPB,  AND 

115  SHOCK  RESULTS 

116  A.  Copper 

117  Figure  2  shows  over  a  larger  range  of  stress  than  in  Fig. 

118  1  a  comparison  of  the  SHPB  results  from  both  Follansbee  et 

119  a/.1  and  the  shock-induced  plasticity  results  of  Swegle  and 


2Uo 
V0  ' 


2kT 

Vo 


In  I  — 


(8) 


148 


in  which  U0  is  the  Gibbs  free  energy  of  activation  in  the  149 
absence  of  an  applied  shear  stress  and  the  constant  e0  is  a  150 
reference  upper-limiting  strain  rate.  The  flow  stress  in  this  151 
case  is  seen  to  be  linear  in  In  e.  152 
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FIG.  3.  HEL  and  shock-induced  plastic  flow  stresses  for  ARMCO  iron  from 
Arnold  (Ref.  4),  in  comparison  with  a  number  of  relationships  proposed  to 
fit  either  slip  or  deformation  twinning  aspects  of  the  strain  rate  dependent 
behaviors. 

^  Equation  (8)  is  plotted  in  Fig.  2,  with  V{)=A*b  taken  as 

154  b3  to  give  2677  V0= 506  MPa  at  300  K.  The  modeled  quan- 

155  tity  2 U0IV0  is  the  theoretical  limiting  stress  for  dislocation 

156  nucleation  and  a  reasonable  order  of  magnitude  estimate  of 

157  —5.8  GPa  was  assumed,  giving  a  value  of  e0  — 2.3 

158  X 109  s-1.  It  is  evident  from  an  inspection  of  Fig.  2  that  the 

159  data  points  of  Follansbee  et  al.  and  of  Swegle  and  Grady 

160  indicate  a  smooth  transition  from  thermally  activated  dislo- 

161  cation  motion  to  thermally  activated  dislocation  generation. 


162  B.  Iron 

163  Much  information  is  available  on  the  shock-induced  de- 

164  formation  properties  of  ARMCO  iron  (first  produced  in  1909 

165  by  the  American  Rolling  Mill  Co.,  ARMCO  ingot  iron  be- 

166  came  the  synonym  for  the  purest  steel-mill-produced  iron, 

167  having  a  purity  of  more  than  99%  iron  and  especially  includ- 

168  ing  very  low  carbon  content)  that  has  been  employed  histori- 

169  cally  also,  as  a  standard  for  equation  of  state  studies,  particu- 

170  larly  relating  to  the  determination  by  Bancroft  et  al.9  of  a 

171  shock  pressure  of  13.2  GPa  for  the  (bcc)  alpha-to-epsilon 

172  (hexagonal  close  packed)  phase  transformation.  Although 

173  previous  shock-induced  deformation  results  were  reported 

174  for  iron  in  several  pioneering  investigations  conducted  at  a 

175  single  grain  size,  Arnold4  investigated  the  shock-induced  de- 

176  formation  properties  for  ARMCO  material  with  the  different 

177  average  grain  diameters  of  20,  40,  80,  and  400  /im,  and,  as 

178  mentioned,  tested  in  plate  impacts  covering  a  large  range  of 

179  projectile  and  target  thicknesses  mostly  centered  on  a  thick- 

180  ness  ratio  of  1/2  for  added  determination  of  the  material  spal- 

181  lation  properties.  In  Fig.  3,  HEF  (closed-triangle)  and  shock 

182  jump  plasticity  (open-circle)  measurements  are  shown  for  the 

183  80  /xm  grain  size  material  that  was  extensively  tested,  for 

184  example,  in  this  case  at  an  intermediate  thickness  ratio  of  3/6 

185  (in  millimeter  thicknesses),  respectively,  for  the  impacting 

186  plate  and  specimen.  Again,  the  shock  pressures  have  been 

187  transformed  to  effective  shear  stress  measurements  according 

188  to  Eq.  (7)  but  with  v=  0.288  for  iron.  Excellent  agreement  is 

189  shown  between  the  open-circle  points  of  Arnold4  and  the  pair 

190  of  filled-circle  measurements  reported  by  Swegle  and 

191  Grady.3 


A  comparison  in  Fig.  3  of  the  power  law  (e1/4)  and  linear 
(In  e)  dependence  given  by  Eq.  (8)  favors  the  latter  descrip-  193 
tion.  Nevertheless,  model  evaluation  in  this  case  of  the  slope  194 
of  the  linear  dependence  on  In  e  in  Eq.  (8),  with  V0 =b3  and  195 
6  =  0.248  nm  for  slip  in  bcc  iron,  gives  the  dot-dashed  line  of  196 
lower  slope  that  is  shown  in  the  figure;  the  slope  is  approxi-  197 
mately  three  times  lower  than  the  fitted  long-dashed  line  198 
through  the  combined  Arnold4  and  Swegle  and  Grady  '  data.  199 
Armstrong  et  al.  accounted  for  the  steeper  slope  of  the  200 
shock-induced  plasticity  results  in  Fig.  3,  corresponding  to  201 
employment  of  V0=(63/ 3),  by  proposing  that  the  propagat-  202 
ing  shock  front  was  nucleating  a  nanoscale  structure  of  de-  203 
formation  twins,  with  V0=b^ji>  and  bT=2(a/6)  [111]  for  a  204 
smallest  two-layer  twin  thickness;10  a  is  the  crystal  unit  cell  205 
lattice  parameter.  Postshock  metallographic  observations  206 
made  on  Arnold’s  80  /x m  grain  size  material  showed  signifi-  207 
cant  deformation  twinning  for  shock  pressures  greater  than  208 
—3  GPa,  corresponding  to  cr>  —  1.8  GPa,  thus  providing  209 
evidence  for  significant  residual  microscale  twinning  associ-  210 
ated  with  all  of  the  open  circle  points.  211 

III.  SLIP/TWINNING  TRANSITION  AT  THE  HUGONIOT  212 
ELASTIC  LIMIT  IN  IRON  213 

An  interesting  comparison  of  predicted  slip  and  defor-  214 
mation  twinning  stresses  also  is  indicated  at  the  lower  level  215 
of  stress  in  Fig.  3  for  the  closed-triangle  cr  values  obtained  216 
from  Arnold’s  HEL  measurements.  In  the  figure,  the  horizon-  217 
tal  double -dot-dashed  line  prediction  was  obtained  for  the  218 
twinning  stress  crT  from  the  equation  219 

(jj-  =  (Tfo  ■+■  ky-f  113  (9)  220 

in  which  crTi)=  330  MPa  and  kT=  90  MPa  mm1/2  as  deter-  221 
mined  from  other  results  reported  for  ARMCO  iron  and  low  222 
carbon  steel  materials,11  whereas  the  Zerilli  and  Armstrong  223 
(ZA)  curve,  that  is  also  shown  in  the  figure,  applies  for  slip-  224 
controlled  deformation  as  computed  with  the  bcc-type  225 
equation7  226 

a  =  B  exp(-  f3T)  +  Ae"  +  (TC  +  ky£~112,  (10)  227 

in  which  /3=/30-/3l  In  e  and  <xc,  B,  /30,  /3h  A,  n,  and  ky  were  228 
taken  from  earlier  reported  experimental  constants  deter-  229 
mined  for  other  results.  The  indication  from  the  compared  230 
relationships  and  the  HEL-transformed  (closed-triangle)  231 
measurements  is  that  the  80  /xm  grain  size  material  may  232 
have  initially  yielded  by  slip  but  only  just  so  because  twin-  233 
ning  occurred  soon  after  deformation  began,  as  definitely  234 
established  via  postshock  metallographic  observations  made  235 
for  the  higher  open-circle  flow  stress  results.4  236 

The  competition  between  slip  and  deformation  twinning  237 
in  iron  and  other  bcc  metals  is  a  well-known  consideration  in  238 
accounting  for  deformation  stress  levels  measured  at  low  239 
temperatures  or  high  strain  rates  near  to  those  required  for  240 
cleavage  fracturing.  Therefore,  it  seemed  worthwhile  to  re-  241 
examine  the  total  HEL  and  flow  stress  measurements  re-  242 
ported  by  Arnold4  for  the  80  /xm  grain  size  material,  as  243 
shown  in  Fig.  4.  The  figure  legend  includes  the  full  range  of  244 
projectile/target  thicknesses  that  were  employed  in  the  tests  245 
and  shows  the  dependencies  of  the  four  constitutive  equa-  246 
tions  that  have  been  described  above;  crSG  is  for  the  Swegle  247 
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FIG.  4.  (Color  online)  A  comprehensive  listing  of  shock-induced  HEL  and 
plastic  flow  stresses  for  ARMCO  iron  at  different  projectile/target  thick¬ 
nesses  from  Arnold  (Ref.  4),  in  comparison  with  a  number  of  relationships 
proposed  to  fit  either  slip  or  deformation  twinning  aspects  of  the  strain  rate 
dependent  behaviors. 

248  and  Grady  power  law  relation  and  crAAZ  is  the  thermally 

249  activated  dislocation  generation  relation,  Eq.  (8). 

250  To  some  extent,  the  distinctions  made  reasonably  clearly 

251  in  Fig.  3  are  now  blurred  by  the  variations  shown  among  the 

252  greater  number  of  experimental  results.  On  the  expanded  ab- 

253  scissa  scale  covering  a  range  of  strain  rates  between  ~6 

254  X  103  and  —4.5  X  106  s-1,  the  lower  HEL  cr  values  appear  to 

255  follow  the  computed  twinning  stress  level  at  first  but  then 

256  move  a  little  bit  upwards  at  the  higher  strain  rates  to  the 

257  calculated  slip  stress  dependence  that  is  indicated  in  Fig.  3. 

258  The  same  linear  in  In  s  dependence  of  the  higher  shock- 

259  induced  plastic  flow  stress  measurements  from  Fig.  3  is 

260  drawn  in  Fig.  4  as  the  long-dashed  line,  and  the  previously 

261  described  data  are  also  shown  to  be  buried  in  the  total  band 

262  of  the  higher  distributed  measurements.  Close  examination 

263  of  the  higher  distributed  stress  values  shows,  at  the  lower 

264  strain  rates,  that  the  shock-induced  flow  stress  jumps  to  a 

265  relatively  constant  plasticity  level  and  persists  so  over  an 

266  approximate  order-of-magnitude  increase  in  strain  rate. 

267  Thereafter,  with  reasonable  variations,  the  shock-induced 

268  plastic  flow  stresses  follow  the  linear  in  In  e  dependence 

269  previously  described  with  Eq.  (8). 

270  IV.  GRAIN  SIZE  DEPENDENT  HUGONIOT  ELASTIC 

271  LIMIT 

272  Despite  the  close  relationship  of  the  slip  and  twinning 

273  stress  measurements  shown  in  Fig.  4  for  the  80  /um  grain 

274  size  material,  the  competition  between  the  two  deformation 

275  mechanisms  is  strongly  influenced  by  the  material  grain  size, 

276  for  example,  as  evidenced  in  the  comparison  of  predictions 

277  from  Eqs.  (9)  and  (10),  and  particularly  involving  the  differ- 

278  ent  magnitudes  of  the  microstructural  stress  intensity,  for  ex- 

279  ample,  of  kT=  90  MPa  mm1/2  for  deformation  twinning"  as 

280  compared  with  a  slip-controlled  yield  stress  value  of  ky 

281  =  22  MPa  mm1/2;  a  lesser  value  of  ke=  5  MPa  mm1/2  applies 

11  12 

282  for  copper.  ’ 

283  Assessment  of  the  relatively  more  important  grain  size 

284  influences  on  the  iron  HEL  cr  values,  therefore,  is  shown  next 

285  on  a  polycrystal  reciprocal  square  root  of  grain  diameter  ba- 

286  sis  in  Fig.  5  for  both  the  slip  ay  and  twinning  oy  stress 
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FIG.  5.  Preshock  hardness  and  shock  flow  stress  measurements,  both  from 
Arnold  (Ref.  4),  along  with  shock  measurements  from  Rohde  (Ref.  13)  and 
Barker  and  Hollenbach  (Ref.  14)  plotted  against  inverse  square  root  of  grain 
size  for  both  slip  and  deformation  twinning. 

dependencies.  Individual  HEL  cr  values  obtained  from  earlier  287 
pioneering  investigations  made  by  Rohde  '  and  by  Barker  288 
and  Hollenbach14  are  included  in  the  figure.  The  importance  289 
of  a  thermally  activated  t*  in  Eqs.  (2)  and  (10)  is  made  clear  290 
by  the  increasing  intercept  friction  stresses  obtained  for  the  291 
c ry  lines  evaluated  at  the  differently  labeled  strain  rates.  For  292 
ay,  a  solid  line  of  correspondingly  higher  kT  was  obtained  293 
from  Eq.  (9),  and  also  a  dashed  line  is  shown  for  a  possibly  294 
raised  ay  at  the  higher  shock-imposed  loading  rates.  The  slip  295 
versus  twinning  issue  for  iron  had  been  assessed  previously  296 
on  this  same  grain  size  basis  for  ARMCO  iron  SHPB  297 
results12  and  for  the  current  material  grain  sizes  subjected  to  298 
Taylor  solid  cylinder  impact  tests.16  299 

Consider  first  in  Fig.  5  the  lower  filled-square  points  that  300 
are  shown  to  follow  an  approximate  ay- type  dependence.  301 
These  points  were  obtained  from  diamond  pyramid  hardness  302 
measurements  made  by  Arnold4  on  the  different  grain  size  303 
materials  before  shocking.  For  the  purpose  of  Fig.  5,  the  304 
individual  hardness  measurements  were  divided  by  a  factor  305 
of  3  to  represent  the  material  flow  stress  at  a  strain  value  of  306 
0.075.  In  turn,  the  lowest  ay  linear  dependence  shown  for  a  307 
conventional  strain  rate  of  0.014  s~  was  raised  by  Acre  308 
=Aen  for  e  =  0.075  by  employing  the  previous  iron  constants  309 
A  and  n  mentioned"  for  Eq.  (10).  Thus,  the  raised  line  pass-  310 
ing  just  above  the  converted  hardness  points  was  obtained.  311 
Reasonable  agreement  is  taken  to  be  shown  by  the  converted  312 
data  points  and  the  slip  comparison.  313 

Next,  attention  is  turned  to  the  higher  filled  diamond  314 
points  for  Arnold’s  converted  HEL  cr  measurements  that  315 
were  obtained  at  different  grain  sizes  but  at  the  same  3/6  316 
projectile/target  ratio.  These  data,  which  are  vertically  spread  317 
over  a  range  in  stress  at  each  grain  size,  are  associated  with  318 
different  strain  rates,  also,  as  indicated  at  the  left-side  two  319 
larger  grain  diameter  cases  of  400  and  80  yum  with  values  of  320 
e—  105  and  106  s-1.  Note  that  the  two  higher  ay  lines  are  at  321 
the  same  pair  of  strain  rates.  Rohde  '  estimated  that  his  filled  322 
circle  point  was  obtained  at  a  strain  rate  equaling  105  s-1.  In  323 
the  Fig.  5  intersection  of  slip  and  twinning  lines  at  the  higher  324 
strain  rate,  then,  the  80  /rm  grain  size  material  is  seen  to  be  325 
very  near  to  the  predicted  transition  from  slip-controlled  326 
plastic  flow  at  smaller  grain  size  and  twinning-controlled  de-  327 
formation  at  larger  grain  size,  particularly,  if  a  higher  twin-  328 
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FIG.  6.  The  strain  rate  dependence  of  the  HEL  yield  stresses  of  ARMCO 
iron  at  different  grain  sizes,  from  Arnold  (Ref.  4),  compared  with  estimated 
twinning  stress  levels  and,  for  a  grain  size  of  80  /xm,  with  the  bcc  ZA 
slip-type  equation  (Ref.  7). 

329  ning  (dashed-line)  stress  dependence  may  be  operable  at 

330  these  greater  strain  rates.  Armstrong  and  Worthington  pre- 

331  viously  reported  a  model  for  deformation  twinning  that  pre- 

332  dieted  a  relatively  small  strain  rate  dependence  of  crT. 

333  There  is  also  indication  in  Fig.  5  of  the  smaller  40  and 

334  20  /xm  grain  size  materials  possibly  following  a  crs  depen- 

335  dence  for  slip.  In  this  way  the  measurements  would  be  fol- 

336  lowing  the  well-known  increasing  tendency  for  slip  to  be  rate 

337  controlling  at  finer  grain  sizes.  Figure  6,  in  showing  details 

338  of  the  strain  rate  dependence,  appears  at  first  sight  to  confirm 

339  the  indication.  In  Fig.  6,  the  abscissa  strain  rate  scale  is 

340  spread  larger  to  show  the  HEL  cr  measurements  for  the  four 

341  grain  size  materials  at  the  respective  plastic  strain  rates  de- 

342  termined  for  the  higher  shock  propagation  stresses.  The  com- 

343  puted  twinning  stresses  from  Eq.  (9)  are  marked  along  the 

344  ordinate  scale.  The  relative  strain  rate  independence  of  the 

345  400  /xm  grain  size  material  and  the  level  of  the  stress  values 

346  compared  to  the  twinning  stress  marker  both  are  in  agree- 

347  ment  with  the  HEL  being  fully  controlled  by  deformation 

348  twinning.  The  dashed  curve  through  the  80  /xm  material  re- 

349  suits  is  that  computed  with  Eq.  (10)  by  employing  the  previ- 

350  ously  determined  constants.  Agreement  is  shown  with  a  slip- 

351  controlled  HEL  cr  but  the  determined  stress  levels,  as 

352  indicated  earlier,  occur  above  the  predicted  twinning  stress 

353  level.  Here  the  indication  is  that  the  80  /xm  material  is  in  the 

354  transition  region  where  slip  and  twinning  may  be  equally 

355  preferred  but  the  measured  cr  values  are  principally  deter- 

356  mined  by  slip.  The  40  and  20  /xm  grain  size  materials  also 

357  appear  to  follow  a  slip-type  strain  rate  dependence  for  the 

358  HEL  cr  but  their  being  higher  than  the  predicted  twinning 

359  stress  level  cannot  be  accounted  for  by  the  slip-controlled 

360  grain  size  dependence  with  ky= 22  MPa  mmI/2,  as  shown  in 

361  Fig.  7. 

362  V.  GRAIN  SIZE  INDEPENDENT  SHOCK-INDUCED 

363  PLASTIC  FLOW  STRESS 

364  Several  models  have  been  proposed  for  dislocation 

365  nucleation  at  a  propagating  shock  front  beginning  from  the 

366  pioneering  description  of  shock  induced  deformations  de- 

367  scribed  by  Smith.  A  main  consideration  in  determining 

368  the  need  for  dislocation  (or  twin)  nucleation  is  that  the  large 
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FIG.  7.  The  HEL  and  shock  stresses  for  ARMCO  iron  of  different  grain 
sizes,  from  Arnold  (Ref.  4),  in  comparison  with  the  differently  modeled  slip, 
deformation  twinning,  and  twin  nucleation  stress  predictions. 

magnitude  of  the  locally  imposed  shear  strains  at  all  points 
along  a  propagating  shock  front  are  too  large  to  be  relieved  370 
by  the  displacement  of  the  few  dislocation  line  segments  371 
contained  in,  or  crossing,  the  front  or  by  the  remote  displace-  372 
ment  of  the  resident  dislocation  density  behind  the  front.  The  373 
intersection  of  grain  boundaries  with  the  shock  front  is  an  374 
even  rarer  occurrence  and  thus,  the  material  grain  size  should  375 
not  be  a  significant  factor  in  determining  the  shock-induced  376 
plastic  deformation  rate.  The  model  consideration  is  in  377 
agreement  with  the  results  shown  in  Fig.  7  for  iron  of  400,  378 
80,  and  20  /xm  grain  size,  all  tested  at  3  to  6  mm  projectile/  379 
target  thicknesses.  380 

In  the  lower  part  of  Fig.  7,  the  just-discussed  HEL  stress  381 
values  for  each  grain  size  are  plotted  in  comparison  with  382 
both  Eqs.  (9)  and  (10)  predictions;  the  40  /xm  grain  size  383 
results  are  omitted  from  the  figure  for  clarity.  As  can  be  seen  384 
at  this  dimensional  scale,  except  for  the  80  /xm  grain  size  385 
material  results,  the  emphasis  would  be  on  a  twinning  expla-  386 
nation  for  the  HEL  stress  and,  as  mentioned  above,  the  dif-  387 
ferent  grain  size  predictions  for  slip  control  are  shown  to  be  388 
too  small  to  account  for  the  results.  The  greater  interest  in  389 
Fig.  7,  however,  is  in  the  apparent  grain  size  independence  of  390 
the  shock-induced  plastic  flow  stresses  at  all  of  the  grain  391 
sizes.  The  double-dot-dashed  line  is  the  same  one  drawn  in  392 
Fig.  3  for  the  twin  nucleation  relationship  given  in  Eq.  (8).  393 
The  result  provides  confirmation  of  the  previously  referenced  394 
shock  model  descriptions.  395 

VI.  ISENTROPIC  COMPRESSION  OF  COPPER  396 

As  mentioned  in  Sec.  I,  the  uniform  pressure  buildup  in  397 
isentropic  compression  provides  for  the  possibility  of  reach-  398 
ing  a  high  plastic  flow  stress  and  corresponding  high  plastic  399 
strain  rate  in  shockless  loading.  Jarmakani  et  al .6  reported  400 
quasi-isentropic  compression  experiments  performed  with  a  401 
two-stage  gas  gun  employing  functionally  graded  impacters  402 
applied  to  [001]-oriented  copper  crystals  sustaining  pressures  403 
ranging  between  17.7  and  51.5  GPa.  404 

Figure  8  shows  the  effective  shear  stress  measurements,  405 
achieved  again  via  transformation  with  Eq.  (7),  for  test  re-  406 
suits  involving  either  zero  hold  time  (short)  (open-squares)  407 
or  10  /is.  (long)  (filled-square)  pressure  pulses.  The  stress  408 
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FIG.  8.  A  comparison  for  copper  of  Hopkinson  bar,  shock,  and  isentropic 
compression  stresses  (Refs.  1,  3,  and  6,  respectively)  as  a  function  of  strain 
rate  and  fitted  with  the  proposed  empirical  and  model  constitutive  relation¬ 
ships  proposed  for  the  respective  strain  rate  regimes. 


values  are  significantly  higher  than  those  in  Fig.  2  but,  per- 

41 0  haps  surprisingly,  the  plastic  strain  rates  are  lower  by  one  to 

411  two  orders  of  magnitude.  For  comparison,  the  previously  de- 

412  scribed  shock  results  for  copper  from  Fig.  2  are  included  in 

413  Fig.  8. 

414  Jarmakani  et  al.6  reported  that  an  amount  of  deformation 

415  twinning,  normally  not  a  major  concern  in  the  high  rate  test- 
41  6  ing  of  copper  except  at  the  highest  loading  rates,  was  de- 

417  tected  in  post-ICE  transmission  electron  microscope  (TEM) 

418  examinations  of  the  specimens  tested  at  the  two  highest  pres- 

419  sures.  As  evident  in  Fig.  8,  the  authors  commented  that  their 

420  results  did  not  follow  the  Swegle-Grady  relationship  for 

421  shocking  of  copper;  rather  the  issue  of  a  slip-twinning  tran- 

422  sition  was  explored  via  a  Preston-Tonks-Wallace  (PTW) 

423  constitutive  equation"  to  compare  estimations  of  critical 

424  pressures  for  twinning  in  both  shock  and  ICE  regimes.  Rem- 

425  tngton  et  a/."  provided  a  comparison  of  similar  constitutive 

426  predictions  made  both  by  PTW  and  ZA  equations  applied  to 

427  the  high  rate  deformation  of  tantalum  material;  see  their  Fig. 

428  1(b).  Also,  it  is  notable  that  the  Jarmakani  et  al.6  TEM  ob- 

429  servations  were  interpreted  to  show  production  of  fewer  dis- 

430  locations  in  the  post-ICE  test  specimens. 

431  In  Fig.  9,  rather  than  consideration  of  a  transition  from 

432  slip  to  twinning,  focus  is  directed  to  the  experimental  linear 

433  dependence  of  cr  on  the  direct  strain  rate.  In  a  previous 

434  report,"4  it  was  proposed  that  the  ICE  and  shock  tests  are 

435  fundamentally  different  in  that  Orowan’s  Eq.  (4)  for  disloca- 

436  tion  generation  is  operative  for  the  shock  case,  as  discussed 

437  above  for  both  copper  and  iron  results,  but  under  the  shock- 

438  less  loading  condition  in  isentropic  compression,  the  strain 

439  rate  is  carried  by  the  movement  of  mobile  dislocations  acti- 

440  vated  within  the  resident  dislocation  density  described  by 

441  Orowan’s  Eq.  (3).  The  additional  complication  in  this  case  is 

442  that  with  a  dislocation  density  of,  say,  ~107  cm-2  and  b 

443  =0.255  nm,  the  average  dislocation  velocity  at  a  shear  strain 

444  rate  of  ~5  X  104  s"1  would  be  —2000  m/s,  compared  to  an 

445  elastic  shear  wave  speed  of  —2900  m/s.  The  estimated  dis- 

446  location  velocity  is  too  large  to  neglect  the  drag  force  on 

447  dislocations  during  their  movement  between  thermal  ob- 

448  stacles  in  the  normal  thermal  activation  description. 

449  Previously,  Zerilli  and  Armstrong  incorporated  disloca- 
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FIG.  9.  ICE  results  for  copper,  after  Jarmakani  et  at.  (Ref.  6),  as  a  function 
of  strain  rate,  showing  linear  stress  dependence  on  strain  rate  in  accordance 
with  model  prediction  for  drag-controlled  slip  (Ref.  24). 


tion  drag  into  a  thermal  activation  model  based  on  Orowan’s 
Eq.  (3)  as  part  of  an  investigation  leading  to  added  drag 
resistance  being  a  discounted  option  for  explaining  the  up¬ 
turn  in  the  SHPB  flow  stresses  shown  in  Fig.  1.  From  the 
previous  investigation,  the  equation  for  a  thermal  component 
of  effective  shear  stress  including  drag  resistance  is 

/  ce  V P'T 

a*  =  B  exp(-  /37j  I  1  -  —  j  I  ,  (11) 

\  Act  / 

where 

c  =  c0m2  fixl  pb2  (12) 

is  a  constant  and  c0  is  the  drag  coefficient  defined  by  the 
equation 

br=c0v,  (13) 

where  br  is  the  force  per  unit  length  of  dislocation  line. 
Equation  (11)  is  an  implicit  equation  for  cr*.  In  the  high 
strain  rate  limit,  the  solution  reduces  to 


(14) 


A  value  of  cl [3\=CQm2 1 pb2~ 0.4  MPa  s  is  shown  in  Fig.  8 
to  fit  the  Jarmakani  et  al.6  data.  The  exact  solution  for  Eq. 
(11)  is  shown  as  the  dash-double-dot  line  in  both  Figs.  8  and 
9.  In  the  expanded  scale  of  Fig.  9,  the  linear  relation,  Eq. 
(14),  is  plotted  as  the  light  solid  line,  which,  in  the  region  of 
strain  rates  of  the  experimental  data,  is  seen  to  be  nearly 
indistinguishable  from  the  exact  solution. 

Determination  of  the  drag  coefficient  from  the  linear 
slope  in  Fig.  9  requires  an  estimation  of  the  dislocation  den¬ 
sity,  which  can  be  made  from  Orowan’s  Eq.  (3).  In  term  of 
the  effective  strain  rate,  Eq.  (3)  becomes 


(15) 
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The  upper  limit  to  the  dislocation  velocity  is  the  shear  wave  478 
speed.  Taking  the  value  of  shear  wave  speed  for  copper  to  be  479 
2900  m/s  and  a  strain  rate  value  e  — 7  X  104  s"1  correspond-  480 
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481  ing  to  the  highest  strain  rate  data  point  of  Jarmakani  et  al.,b 

482  then,  from  Eq.  (15),  a  value  of  p~  2.9  X  107  cm-2  is  deter- 

483  mined.  The  value  is  low  but  not  unreasonable  for  the  mobile 

484  dislocation  density.  With  this  value  of  p,  using  Eq.  (12),  the 

485  estimated  drag  coefficient  cq~  8  X  1 0  4  Pa  s. 

486  Lastly,  there  is  the  issue  of  comparing  the  foregoing  drag 

487  coefficient  to  experimental  values  that  have  been  reported 

488  previously  for  the  movement  of  individual  dislocations,  gen- 

489  erally,  at  orders  of  magnitude  lower  stress  levels.  Jassby  and 

490  Vreeland"fi  reported  such  drag  coefficient  measurements  for 

491  copper  to  be  in  the  range  of  1  X  1CT5  to  8  X  1CT5  Pa  s.  Their 

492  own  measurement  gave  a  value  of  1.7  X  10-5  Pa  s  at  296  K. 

493  The  value  calculated  here  exceeds  the  reported  values  by  at 

494  least  an  order  of  magnitude.  However  it  is  well  known  that 

495  as  the  elastic  shear  wave  velocity  is  approached,  the  disloca- 

496  tion  drag  coefficient  should  increase  without  bound,  becom- 

497  ing  infinite  at  the  shear  wave  speed  (for  example,  see  De 

498  Hosson  et  air1).  Thus,  the  present  result  is  quite  consistent 

499  with  the  idea  that,  at  the  highest  strain  rates  achieved  by 

500  Jarmakani  et  a /./’mobile  dislocation  velocities  are  approach- 

501  ing  the  shear  wave  speed  with  flow  stresses  due  to  drag  ap- 

502  proaching  the  theoretical  yield  stress  limit  for  velocity  con- 

503  trol. 

504  VII.  CONCLUSION 

505  The  strain  rate  dependence  of  the  flow  stress  of  copper 

506  has  been  tracked  over  a  very  large  range  from  ~10-4  to 

507  nearly  107  s-1.  Along  the  way,  a  number  of  different  dislo- 

508  cation  mechanics-based  mechanisms  have  been  associated 

509  with  the  measurements,  for  example,  of  dislocation  intersec- 

510  tions  at  conventional  tension  or  compression  measurements, 

511  of  transition  at  the  high  end  of  SHPB  measurements  to  dis- 

51 2  location  nucleation  that  is  shown  to  be  rate  controlling  for 

51 3  impact-induced  shock  wave  propagation,  and  lastly,  of  ICE- 

514  type  drag-controlled  very  high  dislocation  velocities  being 

515  operative  at  exceptionally  high  stresses. 

516  For  the  higher  ICE  results  displayed  in  Fig.  8,  deviations 

517  from  the  normal  thermal  activation  flow  stress  due  to  drag 

518  begin  appearing  at  strain  rates  as  low  as  500  s_1  and  become 

519  significant  by  1000  s-1.  On  the  other  hand,  the  lower  stress 

520  dislocation  generation  control  shown  for  the  shock  cases  in 

521  Figs.  2-7  tends  to  limit  the  increase  in  flow  stress  to  smaller 

522  values  characteristic  of  an  increased  dislocation  density  and 

523  defers  any  possible  influence  of  drag  to  higher  imposed 

524  strain  rates. 

525  Finally,  the  strain  rate  dependence  of  the  HEL  yield 

526  stress  of  iron  of  different  grain  sizes  tested  in  plate  impact 

527  tests  has  been  accounted  for  in  terms  of  competition  between 

528  alternative  grain-size-dependent  slip  and  twinning  deforma¬ 


tion  responses  but  followed  by  grain-size-independent  828 
shock-induced  plasticity  controlled  by  the  nucleation  of  de-  530 
formation  twins.  The  thermal  activation  areas  determined  531 
from  the  strain  rate  sensitivity  (Act/ A  In  e)T,  as  evaluated  for  532 
twin  nucleation  in  ARMCO  iron  and  for  slip  nucleation  in  533 
shocked  copper,  are  importantly  attributed  to  needed  mecha-  534 
nisms  of  relief  of  the  very  large  strains  imposed  at  all  points  535 
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